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I Abstract. We present some exact solutions and a phase space analysis of metric /(7?)-gravity 

pH , models of the type R". We divide our discussion in « 7^ 2 and n = 2 models. The later model is 

d ' a good approximation, at late times to the f{R) — |:/?tan^' (/?/j3^) gravity model, being this an 

example of a non-singular case. For « 7^ 2 models we have found power law solutions for the scale 
Cf^ ' factor that are attractors and that comply with WMAP 5-years data if n < —2.55 or 1.67 < n < 2. 

I On the other hand, the quadratic model has the de Sitter solution as an attractor, that also complies 

with WMAP 5-years data. 
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INTRODUCTION 



^ . The standard model of cosmology (SMC), which is based upon the Friedmann- 

^ I Robertson-Walker (FRW) models in the context of General Relativity (GR), is suc- 

I cessfuUy to explain the main dynamical and kinematical features of our Universe, such 

as the stretching of light curves of type la supemovae at high redshifts (z ~ 0.3 — 1) 
li, 0, H], CMBR anisotropy measurements [0, and structure formation features 
^ I [la, [70, among other tests. To achieve all this it is necessary to introduce two unknown 

O ■ universe's components: dark matter (DM) and dark energy (DE). It is an unfortunate 

fact that we do not yet have any proved particle physics model that can account for these 
dark components, though some good candidates exist. 
■ In the SMC one employs the Einstein-Hilbert action 



r 1 ^ 4 

Seh = / —RV-gd X, 



(1) 



where R is the Ricci scalar of the space-time, g is the determinant of the metric and 
K = SnG. This gives the well known Einstein field equations 

G^v =R^iv — :j^8^ivR = k{T^v + Tliy^^), (2) 

where T^v includes all the known matter fields and, if necessary, the unknown DM. Posi- 
tive acceleration is possible in virtue of the tensor, which might be a cosmological 
constant or a scalar field, generically called quintessence Js], to which one imposes the 
condition on its pressure to be negative, p < —(1/3) p. The technology known for the 
inflationary dynamics can then be applied to the recent cosmological acceleration [|9|]. 



Alternatively, in recent years it has appeared in the literature a new kind of models, 
generically called /(^, which were initially aimed to account for the present accelera- 
tion of the universe lIlOl lllll. hence avoiding de need of DE. It is possible that the above- 
mentioned observations indicate that GR needs to be modified at low energy scales that 
cosmologically correspond to recent times. More recently these models have been used 
to explain rotation curves and cluster dynamics without DM [14]. 

With the course of the years many models and properties have been published in 
which different fundamental variables has been chosen to vary the lagrangian, i.e., apart 
of the metric, one can also choose the Palatini formalism in which the connections 
are field variables as well, a formalism called Palatini f{R) gravity, in contrast to the 
metric f{R) gravity that only uses the metric. Both of these formalisms use matter fields 
independent of the connections, but one is not forced to do it so. Matter fields that 
possibly depend on the connections are called metric-affine f{R) gravity; for more on 
this topic see [isi 16]. The present paper assumes the standard metric f{R) gravity to 
analyze some exact solutions and to study the phase space. First, we study 7?" gravity 
models with n ^2 and, secondly, we consider a non-singular f{R) model that at late 
times is well approximated by R^. 



f{R) FIELD EQUATIONS 

In order to geometrically modify GR, one can use some of the known gravity invariants, 
and among all possible ones, the simplest modification is to replace the lagrangian ([T]), 



L = R, with a generic function of the scalar curvature [|17LI18|] . f{R)- Thus, the modified 
gravity action is generically written as 

l-f{R)^d\. (3) 

By performing the variation of the action one obtains the field equations of f{R) 
theories [il9i1 : 

/(i?)i?^v-^W(^) +UMvV^Vp-v^Vv)/(i?) = jcr^v. (4) 

This is a coupled nonlinear partial differential equations system of fourth order in the 
metric. The trace of equation dH) is 

3V^V,,/(i?) +Rf'{R)-2f{R) = kT, (5) 

while the trace equation of GR is R = —kT . This is a substantial difference, since 
equation ([5]) is a dynamical law of the scalar curvature, while its counterpart in GR 
is an algebraic law. Thus, we have introduced a new dynamical degree of freedom as we 
have changed the lagrangian from R to f{R) . 

The flat cosmological Friedmann equations for f{R) theories are, using R = 6{H + 
2//2) = 0, 



(6) 



and 



{H + 3//2)/ - 1/ - IHdQf' - dodof = Kp. (7) 
In the following we analyze these equations for some f{R) models. 



R" MODELS 

First, we want to analyze the simplest possible models in the f{R) theories context, 
which are power laws of the scalar curvature, 

/(i?) =Z7-2(«-i)^«^ (8) 

where b is a constant with units of mass and n is a nonzero real number. By using the 
Friedmann equation Q we obtain 



u2{n-\) f.n 

iH + 2H^y-^ {n{n-l)HH + {l-n)H^ + 

{Ar?-ln + A)HH^ + 2{2-n)H^} = Kp. (9) 

First we note that the case R = 6{H + 2H^) = is meaningless for n <2{n I), in 
contrast to GR where i? = in vacuum or for matter fields with T = 0. In vacuum, with 
H + 2H^ ^ 0, equation ^ reduces to 

n{n - l)HH + (1 - n)H^ + {Ar? -ln + A)HH^ + 2(2 - n)H'^ = 0. (10) 

We seek for power law solutions a{t) <x t'^,H = ajt. Using equation (flOl) we obtain 
theequation2(2-n)a2_(4n2_7n + 4)a + 2n2-3n + l =0. One of its roots is a = 1/2 
which is a solution for any n, but it complies with H + 2H^ = and for n < 2 (n 7^ 1) it 
is not acceptable in this description. The other root is given by 

2n^-3n+l 

a = — . (11) 

2 — n 

For an accelerated universe, a > 1, this equation implies 1.36 < n < 2 or n< —0.36. 
The solution (fTT)) is GR equivalent with a dark energy perfect fluid component with an 
effective equation of state parameter given by 

"'//--' + 3(2„'-"lKl-«) - 

The 5-year WMAP data constrains Weff to -0. 1 1 < 1 + Weff < 0. 14 [H]. This implies 
that in a power law model, with n 7^ 2, n is constrained to n < —2.55 or 1.67 < n < 2. 
Note as well that Wgff > ~l can only be achieved for a power law expanding universe. 

By making a Legendre map of the lagrangian density from the nonlinear frame [ 20l ] 
-also called matter frame ifllfl - to the Einstein frame, Carrol et al. showed that a term 



FIGURE 1. Plots V vs x. Here we show results with vq > 1 /2, which asymptotically tend to power law 
solutions that are given by v{x) =const. 



I /If" in the lagrangian can drive an accelerated expansion a{t) <^ t 



(2ra+l)(m+l)/(m+2) 

which is in perfect agreement with equation (fTT)) . where n = —m. 
These power law solutions are attractors, as it can be shown by defining the variables 
x= -H{t),y =H{t) and v{x) = -x^/y^^. Equation ^ then reads 

n(n-l)^—- (2n^ - 3n + 1) ^ + (An^ -ln + A) — + 2(n- 2)/ = (13) 
v-^ ax V 

and an asymptotic behavior to a power law is identified with v{x) ^ a = constant, 
as // — s> 0. By performing numerical integrations for different initial conditions vq we 
obtain the plots shown in Fig. [H for n — —\ and for various values of n that are in the 
interval allowed by the WMAP data mentioned above. Accelerated solutions are found 
for vo > 1/2. Note that vq = 1/2 corresponds to i? = 0. It can be shown that if vq < 1/2, 
the solutions tend to v = 0, i.e. d{t) 0. 



OTHER MODELS 



There exist a plenty of f{R) proposals [16], but have to be analyzed in order to verify 
they pass standard gravitational tests, such as stability criteria or to have a reliable 
Newtonian limit. For instance, theories with f"{R) < and |/"(^)| ^ 1, such as 1 /R"'^ 
or Ini?, are ruled out due to instabilities that appear when trying to model gravity fields 
in stars fl^. 



The initial f{R) models, and some others, add to the lagrangian L = R singular 
i?— functions such as —1/7?'" (with m > 0) having a difficult interpretation to obtain 
a Newtonian limit, since as 7? — > the theory diverges. We think a less restrictive type of 
theories is to consider m < that are non-singular in this sense. 

Motivated by the these facts we propose the following model: 

f{R) = -Rtan-'{R/l5^), (14) 
n 

in which /3 is a constant with units of mass. Indeed, this model passes the stability 
criteria and has reliable Newtonian limit; these matters will be treated in a separate 
work. Fori? ^ [5^ the model behaves basically as f{R) R — f jS^, which is GR with a 
cosmological constant. We define tp as the time when R and we have a transition 

period. To adjust to observations we demand ^ Hq lO^^^eV. 

At late times, R < jS^, the model is well approximated by f{R) ^ 27?^/ n^^. Using the 
Friedmann equation ^ in vacuum we obtain 

2HH + 6HH^-H^ = 0. (15) 
By integrating once, it yields, 

H + 2H^-CH^/^ = 0, (16) 

with C an integration constant. If C = we obtain a t^^-^ and a = const, which are 
not interesting for our aim. In Fig. |2] the phase space is shown for several values of 
C. The long-dashed curve is the equation H = —H^, and the acceleration condition is 
accomplished in the region where the solutions are above it. The double-dashed curve, 
H — 6H^, is the locus of the maximae which occur for all the solutions with C > 0. 
Accelerated solutions occur only for C > 0. The initial condition problem is to specify 

1 /2 

Hq and Hq and then, by using (fT6l) . one determines C = {Hq + 2Hq)/hJ to solve the 
dynamics. 

In Fig. |3] we show the solution with C = 1. In general for C > 0, a phase point A 
with Hq > and Hq < complies with H > 0, as it follows from equation (fT5l) . and 
the solution evolves to the phase point Ha- For a point B with Hq > and Hq > 0, 

> if // > 6//2 and // < if // < 6H^. That is, for all C> 0, the solutions evolve to 
H ^ Ha = C2/^/4i/^, that corresponds to an exponential de Sitter solution. 

We can write equation (fT6l) as 7? = 6CH^/^, and we point out that if i? > the solution 
will have an accelerated final stage. Thus we divide the phase space in three regions: 

I: i? < 0; The solutions are never accelerated and tend to a{t) ^ ao = cte. 

II: The region between R = {H = —2H^) and H = —H^. Here solutions are not 
accelerated but they reach tp in a finite time and enter region III. 

Ill: The solutions are asymptotically de Sitter, a{t) aoe^"' . 

In order to adjust to observations we have ~ //(?|3) ~ //q ~ j8 ~ lO^^^eV, and from 
H{tp) = C2/3 it follows that C ~ p^/^ ~ IQ-^^eV^/^. 

Summarizing, the Arctangent model describes a universe which initially evolves such 
as the one described by GR with a cosmological constant into a stage dominated by 



H 




FIGURE 2. Phase space showing accelerated solutions for curves that are above the long-dashed curve. 



a quadratic term of the scalar curvature. In this stage the solutions tend to a de Sitter 
space-time. 

CONCLUSIONS 

In the phase space analysis presented in this contribution we have analyzed R" and the 
Arctangent model. For the former model with nj^lwe have found power law solutions 
that resulted to be attractors and that comply with the constraints imposed by WMAP 
dil] if n < —2.55 or 1.67 < n < 2, yielding Weff > — 1. The Arctangent model, that was 
motivated as an example of a non-singular /(i?) -gravity, is approximated at late times 
as a quadratic function of R. The resulting equations can be once integrated, leading 
to the simple equation (fT6l) with a constant (C) that depends on the initial conditions 
(Hq, Hq). We have found three phase space regions having the following properties: In 
region I, that corresponds to negative values of C and hence to R < 0, the solutions are 
non-accelerating. In region II, delimited by two parabolic functions, solutions are not 
accelerated but they tend to region III in a finite time of the order of tp , where solutions 
are accelerated and tend to a de Sitter solution. Therefore, regions II and III of the phase 
space end with Wgff = — 1, that is also in accordance with the WMAP results. 
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FIGURE 3. Phase space for C = 1 . Notice that //„ is an attractor to all solutions with C > 0. 
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